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I. 


INTRODUCTION 


The shift towards ever higher frequencies for comiunication links 
and radars has focused attention on the effects of rain in these appli- 
cations. The presence of liquid water in the fonn of rain can have 
deleterious effects in such cases with the major problem being signal 
attenuation which causes deep fading. Hence, a knowledge of the attenu- 
ation to be expected is essential to the design of such systems. 

Weather radar systems at frequencies of high attenuation possess certain 
advantages such as hioh resolution and sensitivity as well as ease of 
calibration using the saturation effect as reported in Reference [1]. 
Again, the measurement of rain rate presupposes a knowledge of the 
reflectivity and attenuation paraiTieters of rain at the frequency of 
measurement. 

Both equivalent reflectivity and specific attenuation can be 
approximated closely as functions of rain rate by power laws of the 
form aR^, where R is the rain rate in nin/hr. Tlieoretical justification 
has been developed in the case of specific attenuation [2] and the 
a and b values have been computed for limited cases by several 
researchers and sunriarized in References [2,3]. However, it appears 
that a consistent tabulation of the a's and b's for both equivalent 
reflectivity and specific attenuation over a broad frequency range does 
not exist in the readily available literature. Consequently this report 
is intended to satisfy this need. 

Tlie results of calculations of equivalent reflectivity, specific 
attenuation and volumetric radar back.catter cross section and of the 
a's and b's for specific attenuation and equivalent reflectivity are 
presented for frequencies from 1 to 500 GHz. Classical Mie theory is 
used to calculate the radar cross sections of spherical drops; and the 
Marshall and Palmer negative exponential drop size distribution is 
assumed. Rain rates from 1.52 to 152.4 nn/hr are used and the tempera- 
ture of the water drops is taken to be 0"C. 


II. CALCULATIONS 


The relevent tl-’oretical results of 
dielectric spheres [4j were recast into a 
computation and are given in tlic Appendix 
the total 
section 


lUon and are given in tnc Appendix, in 
al attenuation cross section, Qt(nm).ai 
I. Qbli^’*")* 4S functions of drop diameter 


Mie for scattering by 
forri more suitable for 
. These are used to calculate 
and the backscalter cross 

D. 


The complex refractive index of water, Oj; , is an important 
parameter in the above calculations . It is strongly dependent on 
frequency and more weakly on temperature. These properties are dis- 
cussed in Reference [5]. 
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Next the Marshall and Palmer drop size distribution was assumed 
[6], It is of the form 

N(D) = e”^^ [drops/m^*m] (1) 

where 

A = aR^ 

N^ = 8000/m^.rrm 
a =4.1 /mm 
e = -0.21 

where R is the rain rate in nri/hr and D is in mm. Then the specific 
attenuation is [7] 


= 4.343 X 10^ j N(D) dD dB/km (2) 

0 


and the volumetric radar backscatter cross section is 


n 


j N(D) dD m^/n? 
0 


(3) 


The equivalent reflectivity may be defined in terms of the 
volumetric radar backscatter cross section as 




6, 3 
mm /m 


where 


K = 


= 10 log^Q 






f lO^A^n 


dBZ 


(4) 


(5) 


( 6 ) 


and A is the wavelength in m. This definition is equivalent to that 
associated with Rayleigh scattering, but remains useful at higher 
frequencies as well. 
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The values of and Zeq were calculated at 36 frequencies from 
1 to 500 GHz for 7 rain rates of 1.52, 2.54, 12.7, 25.4, 50.8, 101.6 
and 152.4 mn/hr. D ranged from .08 to 10.5 mm and the integration 
increment, dO, was .08 nm. The results were then approximated by 
power law equations of the form 

^RD 

= flpQ R dB/km , (7) 


Zgq = a^ R ^ mm^/m^ (8) 

using logarithmic regression analysis to obtain values of a^^Q, b^p, 
a^ and b^. 


III. RESULTS 

Figures 1 and 2 show the equivalent reflectivity, Zeq, and the 
specific attenuation, a^, as functions of rain rate for several fre- 
quencies. The calculated values (represented by crosses) and the power 
law regression lines are shown. Clearly, the deviation from the power 
law is quite small, justifying its use under the assumptions made in 
these calculations. 

The values of ap[), a^ and bop, b^ are shown as functions of 
frequency in Figures 3 and 4 and Table 1. They are seen to be smoothly 
varying functions of frequency. 

Tables 2 and 3 give tl>e calculated values of the equivalent 
reflectivity as a function of frequency and rain rate in mm^/m^ and 
dBZ, respectively. The specific attenuation (dB/km) is similarly 
tabulated in Table 4. Table 5 shows the volumetric radar backscatter 
cross section in m^/m^, also as a function of frequency and rain rate. 

Calculations were also performed from 500 to 1000 GHz. However, 
at the higher frequencies the reflectivity appeared to be quite sensitive 
to the method of numerical integration. This appears to be a result 
of the fact that the incremental drop size, dO, of 0.08 mri could 
be too coarse since it represents a significant fraction of the wave- 
length of 0.3 mn at 1000 GHz. Smaller incremental drop sizes caused 
other computational difficulties; hence, these results were considered 
to be unreliable and were not included in this report. The specific 
attenuation did not exhibit the sensitivity described above for the 
reflectivity. 
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Figure 4. 
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Table 1 

POWER LAW COEFFICIENTS 


PRfO. 


ZEO 


alpha 


IGHZ) 

A 


B 

A 

B 

1 .0 

2.9622C 

2 

1.466 

1.0705E .4 

.605 

1.5 

2.9629C 

2 

1.463 

1.9762E -4 

.666 

?.o 

2.9660E 

2 

1.457 

3.3h94r .4 

,904 

?.5 

2.9724E 

2 

1.450 

5,^467E -4 

.941 

5.0 

2.9717E 

2 

1,442 

7,5l5or ,4 

.963 

5.5 

2 ^ 9336E 

2 

1,444 

1.02llf .3 

1.029 

u.O 

2.eb4eE 

2 

1.454 

1. 34460 .3 

1.071 

^.0 

2.7069E 

2 

1.466 

2.21590 -3 

1.132 

£>,0 

2.t,G22r 

2 

1.323 

3.47o3£ -3 

l.lb3 

7.0 

2.5693E 

2 

1.549 

5,21290 .3 

1.172 

P.O 

2.5936E 

2 

1.564 

7.47110 -3 

1.169 

9.0 

2.fe599E 

2 

1.569 

1.02610 -2 

1 ,lbO 

IP.O 

2.7555E 

2 

1.567 

1,35720 -2 

1.150 

11 .0 

2.6715^ 

2 

1.360 

1.738?r -2 

1,139 

12.0 

5 . 0 0 1 «♦ E 

2 

1.549 

2.16730 .2 

1.13c 

I'.O 

3.4293E 

2 

1.498 

3.73510 .2 

1.106 

20.0 

*♦.096* E 

2 

1.367 

7,36590 -2 

1 .Ceo 

25.0 

4.51.39E 

2 

1.271 

1.23450 >1 

1.046 

iO.O 

‘♦.5969E 

2 

1.164 

1.9496E -1 

i.uio 

i5.0 

3P47E: 

2 

1,070 

2, 6176*- -1 

.972 

HO.O 

5.9775E 

2 

,989 

3.6292E -1 

.934 

50.0 

it.9625E 

2 

.66? 

6,09960 -1 

,6b9 

feO. 3 

2.0507E 

2 

.772 

6,41060 .1 

.619 

70.0 

1.3769E 

2 

.708 

1,05320 0 

.761 

60.0 

9.2273E 

1 

,66? 

1.25710 0 

.752 

90.0 

b . 2254E 

1 

,b29 

1.59160 0 

, 73i> 

10ft. 0 

4,2625E 

1 

.605 

1.52030 0 

.713 

110.0 

2.9707E 

1 

.566 

1.6267f 0 

.699 

l*fC.O 

2.1095E 

1 

.375 

1.71320 0 

.667 

150.0 

6.4153E: 

0 

.556 

1.90460 0 

.663 

?or .0 

2.4546E 

0 

.349 

2.07920 0 

.b69 

230.0 

9,3£43E 

-1 

.351 

2.16730 0 

, 6< 3 

. 500.0 

H.3252E: 

-1 

.554 

2.20940 0 

.617 

350.0 

2.2732E 

-1 

.555 

2,2?4to U 

.fell 

4UC • 0 

1.3312E 

-1 

,553 

2.22230 U 

.bob 

500.0 

5.4676E 

>2 

.550 

2.19660 0 

.605 
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Table 2 



FRro. 

(bHZ) 

1.0 

1.5 
?.0 

2.5 

5.0 

5.5 

4.0 
•S.O 

6.0 

7.0 
A.O 

9.0 
10.0 
11,0 
1?.0 

15.0 

20.0 
^5.0 

30.0 

35.0 

40.0 
bO.O 

60.0 

70.0 

60 . 0 
90.0 

100.0 

110,0 

120,0 

150.0 

200.0 

250.0 

300.0 

550.0 

400 . 0 

500.0 
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Rain RATr 


1. 

27 

2. 

54 

12. 

70 

25 

4.19E 

2 

l.lbE 

5 

1 ,23r 

4 

3.41E 

4. let 

2 

l.lbE 

3 

1.23t 

4 

5.3flE 

4,17E 

2 

1.15E 

3 

i.2ir 

4 

3 34E 

• 

4.15E 

2 

l.lSE 

5 

i.2or 

4 

5,?6E 

4.13E 

2 

I.14t 

5 

1 , IbF 

4 

5,2lf 

4. lit 

2 

l.l3t 

3 

1 , IbE 

4 

3. 15E 

4.09E 

2 

l.l2t 

5 

1,1 jE 

4 

5, lOE 

4,04t 

2 

l.lOE 

3 

1.1 3F 

4 

3.13E 

4.00t 

2 

l.09t 

5 

1.15E 

4 

5,32f 

S.96t 

2 

1.09E 

3 

1.21E 

4 

5.b0r 

3.99t 

2 

1.1 1> 

3 

1.29E 

4 

5.91E 

4.u2t 

2 

1,13E 

3 

1.57F 

4 

4,19L* 

4.09E 

2 

l.lbt 

3 

1.45E 

4 

4 ,42F 

4 ,17E 

2 

1 ,«»0E 

3 


4 

4.bOE 

4.2bt 

2 

l.iJ4t 

3 

l.5oF 

4 

'♦,72E 

4,57t 

2 

1,35E 

3 

1 ,bbr 

4 

4,79F 

5.01E 

2 

1.47E 

3 

l.blF 

4 

4,19t 

5.2uE 

2 

1.47t 

3 

1.57F 

4 

5.25F 

5. ORL 

2 

l,3bt 

.3 

1 .ObE 

4 

2, 'bt 

4.72E 

2 

l.^CE 

3 

B.lwE 

3 

1.65E 

4,21t 

2 

1 ,01E 

3 

5.9bE 

3 

1.15E 

5. loE 

2 

r.7lE 

P 

3.15F 

3 

5.55t 

2.14L 

2 

4,2bt 

2 

1 ,b9F 

3 

2.61E 

1.45t 

2 

2.7lt 

2 

9.44F 

2 

1.50E 

9,77» 

1 

1 , 74f 

2 

5.5oE 

2 


b.b5E 

1 

1.14E 

2 

3. 3bF 

2 

5.07F 

4. bet 

1 

7.bOt 

1 

2.13E 

2 

5,l7t 

5,23E 

1 

5.21E 

1 

1 ,4uF 

2 

^.07t 

2,31t 

1 

3,b5t 

1 

9,53F 

1 

l,40E 

9, 3eE 

0 

1.43t 

1 

3.53F 

1 

5, 14F 

2.7RE 

0 

4 ,12E 

0 

9,93f 

0 

1 ,45E 

l.U7t 

0 

l.5bE 

0 

3, 7bF 

0 

5,53t 

4.9HL 

.1 

7.27t 

-1 

l,7bF 

0 

2.5bF 

2,b3t 

.1 

3.83E 

-1 

9, 19F 

.1 

1.35E 

1.54E 

.1 

2,*:4r 

-1 

5.53F 

• 1 

7.65F 

b . 3bt 

.2 

9, 15r 

-2 

2. IbE 

• 1 

3,17c^ 


40 50.60 101.60 152.40 

4 9.42E 4 2.59E 5 4.66E 5 

4 9.31E 4 2.55L 5 4,56E 5 

u 9 i**r 4 I 4«»t. 5 4 ^h5L 5 

4 6*92E 4 2*4(jE 5 5 

4 eIfe7E 4 2l32t 5 4,l0t 5 

4 b,49f 4 2.30E 5 4,l5t 5 

4 6,4far 4 2,37t 5 4.40L 5 

4 6.9bE 4 2.66E 5 5.16L 5 

4 9,94f 4 3, OPE 5 5,97t 5 

4 l.lOe 5 3.42t 5 b.58L 5 

4 1.20E 5 3.67E 5 b,9fet 5 

4 1.28E S 3.82L 5 7,14t 5 

4 1,33^ 5 3,69t 5 7,15 l 5 

4 l.IbE 5 3,eeE 5 7.02E 5 

4 1.3TE 5 3.61E 5 6.7bL 5 

4 1,50F 5 3.37E 5 5.70E 5 

4 1,o3E 5 2.35E 5 3.70E 5 

4 7.21E 4 1.50E 5 2,23t 5 

4 4,82E 4 9.24E 4 1.32L 5 

4 3,lfaE 4 5.70E 4 7,84t 4 

4 2. ORE 4 3.56E 4 4,78E 4 

3 9,?9r 3 1.49E 4 1.93L 4 

3 4,4bE 3 b.e7L 3 6,72 l 

3 2. Sir 3 3.4bE 3 4.35L 

2 l.2«f 3 1.69t 3 2.3bt 

2 7.53F 2 l.lOE 3 1,37E 

2 4,b7E 2 b,79L 2 e,H4t 2 

2 3,o 3E 2 4,39E 2 5.4bE 2 

2 2.09E 2 2,96E 2 3,b7L 2 

1 7,47r j 1.09E ** 1 . 35E 2 

1 2. HE 1 3.09E I 3.ee£ 1 

0 b.l2E 0 1.20E I l.SlL 1 

0 3.7»E 0 b.blE 0 7.09L 0 

0 2.00E 0 2.97E 0 3.7bE 0 

-I l.lbr 0 1.73E 0 2.l9t 0 

-1 4,b9f .1 7.00E -I 6.tJ9E -I 
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IV. SUrtlARY 


The equivalent reflectivity, Zeq, the specific attenuation, 
and the volumetric radar backscatter cross section, n, associated 
with rain rates from 1.27 to 152.4 m/hr were calculated for 36 fre- 
quencies from 1 to 500 GHz. The rain drops were considered to be 
dielectric spheres and classical Mie theory was used. No corrections 
for distortion of the drops were made. The Marshall and Palmer drop 
size distribution was assumed and tne rain ten’ ,‘rature was taken to be 
0®C. 


The results for Zeq, and were then shown to be closely approxi- 
mated by power laws of the form aR^ over the entire frequency and rain 
rate ranges considered. The values of the a's and b's were given and 
shown to be smoothly varying, well behaved functions of frequency. 
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APPENDIX 


SCATTERING BY A DIELECTRIC SPHERE 


Consider a dielectric sphere having a radius a and complex 
relative dielectric constant 



where , 

= relative complex dielectric constant 
c = real part of c 

r 

c = dielectric constant of free space 
0 

0 = electrical conductivity of the dielectric 

u) = radian frequency. 


(Al) 


The complex index of jefraction is 


n = ► e 
c r 


Let a plane wave 


V = E a 

o X 


(A2) 


(A3) 


be incident on this sphere. 

= incident electric field of wavelength X 
k ~ ^ - propagation constant in free space. 

Then the attenuation cross section, Q^, is [4] 

r “ 1 

I (2"*'>(a + b )] (A4) 

k J 

and the back scatter cross section, is 
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k ‘n=l 


(Ab) 


where the Mie coefficients and b are niven by 

n n ' 


jp(n^p)[pjj^(p)]'-j^(p)[ri^pj^(n^p)] ' 

"" ’ l\'/'(p)[n,pj„(n^p)]'-,i^^(n^p)[pl.<^'(p)] 


(A6) 


b = j > ["c** jp <"c^J ' '”c jn ‘y ) [a j „ (P ) ] ' 

" n^jn(ricP ) [ph^J^ > (p )] ' ' (p )[n^p (n_.p ) ] 


(A7) 


Jn spherical Hankel and Bessel and Neumann functions, 

respectively, where 


P = ka 

= j„(p) - .1 M^(P) 


(A8) 

(A9) 


Lettinq or h^^ and usinq the recursion relationships [7] 


dp ° H+T 


(AlO) 


and 


d[pZ„(p)] 

3p 


p z,;(p) 


* Zn(^) = 2„(P) 


2n+l 


n-1 


(p). 


the derivatives can be eliminated from Lquations (A6) and (A7). For 
nM a^^ and b^ can be rewritten in the fonn 


lb 





- n^j^(p)[nj^_,(n^p) - (n+1 )j^^, (n^p)] , 

r ^ 

Khr(p)[n,i^_,(n^p) - (n.l)j^^,(n_.p)] 

- j„(n^p)[nh<^j(p) - (n+l)h<2)(p)] 



(.M2) 



which is more suitable for computational purposes. 
^ / 


Similarly, for n=0: 


(A14) 
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cjo(ncP)[nc*'o^’(p) - 

t'*^*(p)[jg(n|.p) - n^pj^(n^p)]\ 


(A15) 
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